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ABSTRACT 

We present a sample of 42 high-mass low-metahicity outliers from the mass-metallicity relation of 
star-forming galaxies. These galaxies have stellar masses that span log(M^/Mo) ~ 9.4 to 11.1 and are 
offset from the mass-metallicity relation by —0.3 to — 0.85dex in 12 -I- log(0/H). In general, they are 
extremely blue, have high star formation rates for their masses, and are morphologically disturbed. 
Tidal interactions are expected to induce large-scale gas inflow to the galaxies' central regions, and 
we find that these galaxies' gas-phase oxygen abundances are consistent with large quantities of 
low-metallicity gas from large galactocentric radii diluting the central metal-rich gas. We conclude 
with implications for deducing gas-phase metallicities of individual galaxies based solely on their 
luminosities, specifically in the case of long gamma-ray burst host galaxies. 
Subject headings: galaxies: abundances - galaxies: evolution 



1. INTRODUCTION 

Star-forming galaxies fall on a luminosity-metallicity 
relation such that more luminous galaxies tend to have 
higher gas-phase metallicities (the proxy for which is typ- 
ically the oxygen abundance in units of 12 -f log[0/II]). 
This relation is observed to hold — shifted to lower 
metallicities — at redshifts as high as z ~ 3.5 (Erb et al. 
2006; Maiolino et al. 2008), and its scatter decreases 
to ~ 0.15 dex at z ~ when galaxy stellar mass re- 
places luminosity in the relation (Tremonti et al. 2004). 
Most commonly-accepted theories as to the origin of 
the mass-metallicity relation have as a central propo- 
sition that low-mass galaxies are metal- deficient, rather 
than that high-mass galaxies are metal-enhanced (Larson 
1974; Dalcanton 2007; Finlator & Dave 2008). In fact, 
there is some evidence that the mass-metallicity relation 
may flatten at large stellar masses, log(Mi,/M0) > 10.5 
(e.g., Tremonti et al. 2004), though it is unclear to what 
extent this flattening is due to a saturation of the metal- 
licity indicator used at high oxygen abundances (see e.g., 
Kewley & Dopita 2002; Bresohn 2006; Kewley & Elhson 
2008). 

Because most star-forming galaxies do lie on a mass- 
metallicity locus, we can also learn about the gas-phase 
metallicity evolution of galaxies by studying the proper- 
ties of galaxies that do not fall on the relation. In Peeples 
et al. (2008, hereafter Paper I), we explored the popu- 
lation of low-mass high-metallicity outliers, and postu- 
lated that these metal-rich dwarf galaxies have low gas 
fractions and are therefore nearing the end of substantial 
epochs of star formation. Here we investigate the other 
corner of the mass-metallicity plane by asking what we 
can learn about the evolution of massive galaxies from 
the properties of the high-mass low-metallicity outliers. 

As shown in Figure 1, we find 42 low-metallicity high- 
mass galaxies with masses ranging from log(M^/M0) ~ 
9.4 to 11.1 and offsets from the central mass-metallicity 
relation of —0.3 to —0.85 dex. We describe in § 2 how 
we selected this sample and verified the galaxies' outlier 



status. In § 3, we describe the physical properties of these 
galaxies and discuss possible origins for their low oxygen 
abundances. Specifically, we find that they have highly 
disturbed morphologies strongly suggestive of merging 
or post-merging systems, are extremely blue, and have 
high specific star formation rates. We summarize our 
conclusions and state some implications of these findings 
in §4. 

2. SELECTING HIGH-MASS LOW-METALLICITY 
GALAXIES 

We began with a sample of ~ 110000 star- forming 
galaxies from the Sloan Digital Sky Survey (SDSS) Data 
Release 4 (Adelman-McCarthy et al. 2006) with gas- 
phase abundances measured by Tremonti et al. (2004) 
and stellar masses derived using the techniques of Kauff- 
mann et al. (2003). High- luminosity outliers from 
the mass-metallicity relation can be found in the low- 
metallicity region of parameter space for one of two rea- 
sons: either they have spuriously low derived metallic- 
ities or they are genuine outliers. (The third possibil- 
ity of an object having a spuriously high measured lu- 
minosity is much less likely.) In general, the strengths 
of the visible-wavelength oxygen forbidden [O II] and 
[O III] emission lines increase as the oxygen abundance 
decreases; hence, it is possible for a bright galaxy with 
a weak active galactic nucleus (AGN) component to ap- 
pear to have a low oxygen abundance due to contribu- 
tion from the AGN itself. Likewise, an otherwise "red 
and dead" galaxy with strong [O II] emission might be 
labeled as star-forming and, in the Bayesian analysis of 
Tremonti et al., be assigned an artificially low metallic- 
ity. We therefore selected our sample of massive low- 
metallicity galaxies by first ensuring that they are out- 
liers in the luminosity- and mass-metallicity parameter 
spaces and then following up with several line-ratio diag- 
nostic tests in an attempt to guarantee that their oxygen 
abundances relative to the rest of the sample are indeed 
believable. Images of the galaxies in our final sample are 
shown in Figure 2 and summary information is provided 
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Fig. 1. — Plot of 12 + log(0/H) v. log(M*/M0). The small grey points represent galaxies from the main Tremonti et al. (2004) SDSS 
sample described in § 2, and the step-curve is the median of these points in bins of 0. 1 dex of total stellar mass. The high-mass low-metallicity 
massive galaxies in our sample are plotted as the large grey circles in the lower-right region of the diagram; the only two clearly spiral 
galaxies in the sample are denoted with "§" symbols. For reference we also plot the high-metallicity outliers from Paper I as large open 

points (upper left). 

in Table 1. As only two of the galaxies in Figure 2 arc 
clearly undisturbed spirals, we note these two galaxies 
in Table 1, Figure 1, and subsequent figures with "§" 
symbols. 



a{u) < 0.2 for M > 19.5, 



(3) 



2.1. Sample Selection 

We summarize the selection process for our massive 

metal-poor galaxy sample in Table 2. We begin by defin- 
ing the parent galaxy sample from which our galaxies are 
outliers. Starting with the entire Tremonti et al. (2004) 
sample of star-forming galaxies, we choose only those 
galaxies with SDSS magnitude errors in the g, r, i, and 
z-bands of < 0.1 mag. In the M-band we take those galax- 
ies that satisfy the conditions 

cr(u) < O.lforu < 18, (1) 
a{u) < (j^^ u + 16-5 for 18 < u < 19.5, and (2) 



where u is the observed w-band magnitude and is the 
SDSS error on u. These cuts yield a parent population of 
86754 galaxies. To find the low-metallicity outliers from 
the luminosity-metallicity locus, we take bins of absolute 
g-band magnitude Mg of width 0.4 mag, and each bin we 
take the 1% of objects with the smallest 12 + log(0/H), 
as demonstrated in Figure 3. Because we do not want to 
preferentially select blue objects, we similarly take the 
1% of galaxies with the smallest 12 -|- log(0/H) in bins 
of Mz of width 0.4 mag. Finally, we ensure that these 
galaxies are also outliers in the mass-mctallicity plane 
by taking the 2.5% smallest 12 + log(0/H) galaxies in 
bins of stellar mass of width AlogM^, = 0.1 dex and the 
2.5% largest logM^, galaxies in bins of 12 -|- log(0/H) of 
width 0.1 dex. These selection criteria yield a sample of 
113 candidates for detailed evaluation and further down- 
selection of the true low-metallicity outliers. 
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Fig. 2. — Normalized SDSS g-band images of massive low-metallicity galaxies scaled to 40 x 40kpc, with mass decreasing to the right 
and down from logM^ fa 11.1 to 9.4. The only two clearly spiral galaxies are the middle two on the top row. 



One possible contamination source of these seemingly 
bright metal-poor galaxies is a low- level AGN, such as 
a low-ionization nuclear emission-line region (LINER). 
While it is possible for galaxies with weak AGN activity 
to also have low oxygen abundances, we choose to not 
try to disentangle these two effects on the emission line 
fluxes. Using the line fluxes measured by Tremonti et al. 
(2004)^, we put these galaxies on the standard Baldwin, 
Phillips, & Terlevich (BPT) 1981 diagrams. In partic- 
ular, as shown in Figure 4, we find that about half of 

^ See http: //www. mpa-garching.mpg. de/SDSS/DR4/raw_data. html 
for this very rich dataset. 



these outliers have strong enough [O I] A 6300 to place 
them in the "AGN" region of the BPT diagram. (We 
note that while the Tremonti ct al. (2004) galaxies were 
chosen to be star-forming based on where they fall in the 
[O III]/H/3 V. [N II] /Ha plane, 3.4% of the parent sam- 
ple galaxies fall in the AGN region of the [O III]/H/3 v. 
[O I] /Ha plane.) As expected, the strong- [O I] galaxies 
also fall close to the "composite"' H II-AGN boundary 
as defined by Kewley et al. (2006b) in the [O III]/H/3 v. 
[N II] /Ha plane. Visual inspection of the spectra and 
images of these bogus outliers reveals that many of them 
are quite red, with only strong [O II] A 3727 and perhaps 
weak Ha emission being suggestive of ongoing star for- 
mation. Because [O I] A 6300 is a relatively weak line and 
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LOW-METALLICITY MASS-METALLICITY OUTLIERS 
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Note. — Sample of massive low-metallicity galaxies, reverse sorted by 
stellar mass. RA and Dec are in degrees (J2000.0), 12 + log(0/H) and 
stellar masses are from Tremonti et al. (2004). 



^Spiral galaxy. 



TABLE 2 
Cuts for Sample Selection 



Cut 



Number Surviving 



Low magnitude errors 86754 

1% low O/H w.r.t. Mg 855 

1% low O/H w.r.t. Mz 598 

2.5% low O/H w.r.t. logM* 381 

2.5% high logM* w.r.t. O/H 113 

within 2(7 of H II region of [O I] BPT diagram 56 

OK spectrum reduction 52 

reproducible [N II] /Ha ratio 44 

2.5% low PP04 (O/H) w.r.t. logM,^ 42 

Note. — See § 2 and Figure 3 for a more detailed explanation. 




Fig. 3. — An example of one of the cuts made to find low- 
mctallicity high-luminosity outliers from the mass-mctallicity re- 
lation. Plotted is 12 + log(0/H) vs. Mg for the main sample {small 
light grey points), the 1% objects the lowest log(0/H) in bins of 
Mg of width 0.4 mag {small black points), and the final sample of 
high-mass low-metallicity outliers {large grey points). The only two 
clearly spiral galaxies in the sample are denoted with "§" symbols. 
The histogram shows the 1% threshold in 12 + log(0/H) in bins of 

Mg. 

many of the outliers arc clustered near the the Kcwlcy 
et al. (2006b) boundary, we kept galaxies whose [O I] /Ha 
ratio is within 2-cr of the H II-AGN boundary and whose 
spectra passed a visual inspection test for, e.g., clear H/3 
emission. Representative spectra are shown in Figure 5. 

We further tested whether or not the measured metal- 
licities could be due to measurement errors in the rele- 
vant line fluxes. Using the Data Release 6 SDSS spectra 
(Adclman-McCarthy et al. 2008), we subtracted the un- 
derlying stellar continuum using the STARLIGHT pro- 
gram (Cid Fcrnandcs et al. 2005), which simultaneously 
corrects for reddening and extinction using the Cardelli 
et al. (1989) extinction curves. While the extinction cor- 
rections are all small {Ay < 1), the Tremonti et al. line 
fluxes do not account for reddening, and thus a compar- 
ison of widely spaced line ratios produces a systematic 
offset. We therefore reject seven galaxies for which our 
measurements strongly disagree with Tremonti et al.'s on 
the [N II] /Ha ratio; the difference in each of these cases 
is such that using the Tremonti et al. ratio would result 
in a lower 12 + log(0/H) than if ours were used. The 
added advantage to using the [N II] /Ha ratio is that this 
is the ratio used in the Pettini & Pagel (2004) metallicity 
determination, which we use as discussed below. 

2.2. Metallicity Believability 

We note that we do not impose a metallicity error cut 
when selecting these low-metallicity outliers. This is be- 
cause the errors on Tremonti et al.'s 12 + log(0/H) are 
not a simple function of metallicity in this range. It has 
long been recognized that accurately measuring low oxy- 
gen abundances (e.g., using R23) can be tricky due to 
degeneracies in the ionization parameter, and the line 
ratios usually used to determine the ionization parame- 
ter (e.g., [O'lII]A5007/[O II]A3727) are also dependent 
on metaUicity (McGaugh 1991; Kewley & Dopita 2002). 
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Fig. 4.— Baldwin, Phillips, & Terlevich diagrams of log([0 III]/H^) v. log([N II]/Ho) [left] and log([0 III]/H/3) v. log([0 I]/Ha) [right]. 
Plotted are the spectroscopic galaxies from SDSS DR4 (small grey points), the galaxies passing our luminosity, mass, and log(0/H) cuts 
but not in the final sample {open points), the galaxies passing the BPT diagram cuts {light grey points), and the galaxies in our final 
low-metallicity sample {dark grey points). Lines are taken from Kewley et al. (2006b). 



Unfortunately, this means that many mctallicity indica- 
tors are discontinuous or systematically tend to avoid 
particular abundance solutions. 

Wc therefore adapted a different approach to ensure 
that our low-metallicity high-mass outliers geniunely 
have relatively low oxygen abundances. As discussed 
in §2, we first chose only those galaxies that lie in the 
H 11 region locus of the BPT diagrams, and secondly, 
we discarded those galaxies for which we did not mea- 
sure a comparable [N II] /Ha ratio to those given by 
Trcmonti ct al.. Finally, we argue that if those galax- 
ies are truly high-mass low-metallicity outliers, then they 
should be outliers regardless of how 12-|-log(0/H) is mea- 
siucd. (See Kcwlcy & Dopita 2002 and Kewley k Ellison 
2008 for thorough discussions of how and why diiferent 
metallicity callibrations differ.) We therefore used the 
Tremonti et al. measurements recalculate 12 -|- log(0/H) 
for the entire parent sample using the Pettini & Pagel 
(2004) method, where the oxygen abundance is given by 

12-t-log(0/H) = 9.37-^2.03 xN2-t- 1.26 xN2^-t-0.32xN2^ 

(4) 

where N2 = log([N II]A6584/Hq!). The main disadvan- 
tage of this method is that it uses nitrogen as a proxy 
for oxygen; however, there are several factors that make 



it advantageous in these objects. First, it is reddening 
insensitive because the [N II] A 6584 line and Ha are sep- 
arated by only ~ 20 A. This means that we can safely 
use the Trcmonti ct al. measurements and our own mea- 
surements without being overly concerned about the red- 
dening corrections. Second, unlike the various of R23 
methods, the Pettini & Pagel method is both continuous 
and single- valued throughout the metallicity range of in- 
terest; there are no "upper-" and "lower" -branches to 
concern us. We show the mass -metallicity relation with 
12-hlog(0/H) measured with the Pettini & Pagel (2004) 
method in Figure 6. Only one galaxy in the low-[0 I] 
sample did not pass the 12 + log(0/H)-logM* cuts de- 
scribed in Table 2, suggesting that the remaining galax;ies 
are true low-metallicity outliers. 

We have made several other checks for potential sys- 
tematics. First, spiral (and thus many star-forming) 
galaxies are known to have metallicity gradients such 
that their centers have higher abundances than at larger 
radii (e.g., Zaritsky et al. 1994; Kcnnicutt et al. 2003). 
It is therefore reasonable to expect that, on average, spi- 
ral galaxies with larger fractions of their light falling in 
the 3" SDSS fiber will have lower integrated metallicities. 
However, most of our galaxies (see Fig. 2) are not spiral 
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Sample spectra. Two low-metallicity high-mass outliers {top and middle panels) and one galaxy rejected as a potential AGN 



galaxies, so we do not expect for this to be a huge effect. 
Regardless, as shown in Figure 7, while the galaxies in 
our sample do have somewhat high fiber fractions, this 
is only a small effect, and the offsets wo sec in mctal- 
licity are much greater than can be explained by large 
fiber fractions, which would cause contamination of cen- 
tral high- met allicity H II regions by low-metallicity H II 
gas at larger galactocentric radii. 

Similarly, the final rcdshift distribution for the galax- 
ies in our sample peaks at higher z than for the main 
sample; since we have convinced ourselves that this is 
not purely a fiber fraction efi^ect, it is worrisome that it 
might partially be a redshift-evolution effect. When the 
universe was younger, galaxies arc expected to have had 
lower oxygen abundances; the observed mass-metallicity 
relations for redshifts z > 0.1 are systematically shifted 
to lower metallicity from what is observed for nearby 
galaxies (Erb et al. 2006; Maiolino et al. 2008). We plot 
our low-metallicity outliers against the observed mass- 
metallicity evolution from Maiolino et al. (2008) in Fig- 
ure 8. While the median redshift for galaxies in our sam- 
ple is z ~ 0.17, the observed oxygen abundances are more 
typical of redshifts 2.2 < z < 3.5. Though a large amount 
of scatter in the mass-metallicity relation is both ex- 
pected and observed (Savaglio et al. 2005; Erb et al. 2006; 




log (M./Mg) 

Fig. 6.— Plot of Pettini & Pagel (2004) 12 -|- log(0/H) v. 
log(M*/M0). The galaxies in our sample are plotted as the large 
grey circles in the lower-right region of the mass-metallicity di- 
agram; the only two clearly spiral galaxies in the sample are 
denoted with "§" symbols. The histogram denotes the median 
12 -|- log(0/H) in bins of stellar mass, while the curves are fits to 
standard deviation offsets as labeled. 



Massive Metal-Poor Galaxies 



7 




0,4 

fiber fraction 

Fig. 7. — Plot of residual in 12 + log(0/H) from the Tremonti 
et al. fit to the mass-metallicity relation v. fiber fraction. The 
high-mass low-metallicity outliers are shown as the large dark grey 
points at the bottom of the diagram, with the only two distinctly 
spiral galaxies in the sample as "§" symbols. For reference we 
plot the high-metallicity outliers from Paper I as large open points 
(top). 
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Fig. 9. — Mass— metallicity outliers on the color-magnitude dia- 
gram [{g — r) vs. Mr, both fc-corrected to z = 0.1 using the tem- 
plates from Asset et al. (2008)] from the SDSS DR6 spectroscopic 
galaxy sample (Adelman-McCarthy et al. 2008). The blue circles 
denote the 42 low-metallicity high-mass galaxies in our sample; for 
reference, metal-rich dwarf galaxies from Paper I are denoted by 
the grey circles. 
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Fig. 8. — Evolution of the mass-metallicity relation from 
Maiolino et al. (2008) with comparison to low-metallicity high- 
mass outliers (grey points). Fits to the mass-metallicity relation 
from Maiolino et al. (2008) are shown with lines: z ~ 0.7, dot- 
ted; z ~ 0.7, short dashed; z ^ 2.2, long dashed; and z ^ 3.5 
solid. The grey points denote our high-mass low-metallicity out- 
liers (0.035 < 2 < 0.27), with re-calculated metallicities using the 
Kewley & Dopita (2002) diagnostic using our measurements of the 
line ratios so that they are on the same scale as the Maiolino et al. 
results. (The stellar masses have also been shifted to match the 
Maiolino et al. fits.) 

Kobayashi et al. 2007; Finlator & Dave 2008; Maiolino 
et al. 2008) at all redshifts, the abundance offsets of our 
low-metallicity outliers are much more pronounced than 
the relatively subtle expectations from metallicity evolu- 
tion with cosmic time. 

3. DISCUSSION 



To seek to explain these galaxies' low oxygen abun- 
dances, we need to consider the population's other phys- 
ical properties. As shown in Figure 9, the galaxies in 
our sample are all very blue; most — but not all — are 
also outliers in the galaxy color-magnitude diagram. In 
§3.1, we discuss how this blueness can be attributed to 
the galaxies' high specific star formation rates. Fig- 
ure 2 shows that 40 out of 42 of the galaxies in our 
sample have disturbed morphologies suggestive of merg- 
ing systems. Though some of the less-resolved galax- 
ies in our sample may appear only slightly morpholog- 
ically disturbed in the SDSS images, evidence suggests 
that, with higher resolution, these galaxies do in fact 
have unusual morphologies. In a Hubble Space Tele- 
scope study of local Lyman break galaxy (LBG) analogs, 
Overzier et al. (2008) found SDSSJ102613.97-h4884458.9 
{[a,S] = [156.5580,48.7497], leftmost image on the bot- 
tom row of Figure 2) to have a strongly asymmetric pro- 
file with several distinct knots of star formation. In a pre- 
liminary study of the most UV luminous galaxies, Hoopes 
et al. (2007) found that these high specific star formation 
rate LBG analogs are metal-poor by ^ 0.5 dex relative to 
other galaxies of similar stellar mass. Hoopes et al. and 
Overzier et al. suggest that the observed high specific 
star formation rates and low metallicities for these LBG 
analogs are related to the galaxy collisions that formed 
these objects. (See also Struck et al. 2008 for a discus- 
sion of so-called "delayed" galaxies, i.e., high star for- 
mation rate interacting galaxies which have apparently 
managed to retain most of their gas until z = 0.) In § 3.2, 
we show how a simple merger-induced gas infiow picture 
can account for an offset in metallicity, and we discuss 
possible shortcomings in this interpretation. The other 
two galaxies in the sample are undisturbed spiral galax- 
ies, and while these are clearly very interesting objects, 
we have no explanation to offer for their extreme offsets 
from the mass-metallicity relation. 
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Fig. 10. — Specific star formation rate vs. stellar mass. The grey 
points are galaxies from the main Tremonti et al. (2004) SDSS 
sample described in § 2, the low-metallicity massive galaxies in our 
sample are denoted with blue circles (with the only two clearly 
spiral galaxies marked by "§" symbols), and the metal-rich dwarf 
galaxies from Paper I are shown with large grey points for reference. 
The histogram denotes the median specific star formation rate as a 
function of stellar mass and the dotted lines show lines of constant 
star formation rate. The star formation rates estimates are from 
(Brinchmann et al. 2004) and the stellar mass estimates are from 
Tremonti et al. (2004) and Kauffmann et al. (2003); both have been 
aperture corrected. 

3.1. Star Formation Rates 

Ellison et al. (2008a) have shown that galaxies with 
higher specific star formation rates at a given mass 
have preferentially lower metallicities. In Figure 10, we 
plot the specific star formation rate (SSFR) as a func- 
tion of stellar mass for the SDSS star-forming galax- 
ies. Most (but not all) of our low-metallicity outliers 
clearly have higher specific star formation rates than typ- 
ical of the larger sample. We plot the specific star for- 
mation rates against the mass-metallicity residual from 
Tremonti et al. (2004) in Figure 11; generally speak- 
ing, aperture corrections for the specific star formation 
rate only matter at the low SSFR end. While our low- 
metallicity outliers do have higher-than- normal SSFR.s, 
this relative difference alone is not enough to explain 
their low oxygen abundances. 

3.2. Observed Metallicity of Merging Galaxies 

As can be seen in Figure 2, most of galaxies in this sam- 
ple are tidally interacting. Simulations have shown that 
galaxy major mergers are accompanied by starbursts and 
gas infiow from large galactocentric radii (e.g., Barnes 
& Hernquist 1992; Mihos & Hernquist 1996; Cox et al. 
2006). Because typical spiral galaxies have metallicity 
gradients (Kennicutt et al. 2003), the infiowing gas will 
have a lower oxygen abundance than the native central 
gas. This suggests that perhaps our galaxies' observed 
low metallicities are a consequence of large quantities 
of low-metallicity gas from large galactocentric radii are 
being transported into the central few kiloparsecs and 
diluting the metal-rich central gas. 

We can test this idea at the order of magnitude level 
to verify that it can give the observed oxygen abundance 
offsets of -0.3 to -0.8 dex (see, e.g.. Fig. 7). The SDSS 



data measure the galaxies' central oxygen abundances; 
we wish to know what the change in oxygen abundance 
would be if low-metallicity gas from larger radii were 
dumped into the central regions. If we assume that all of 
the gas out to some large radius R was to participate in 
this inflow, then the lower abundance would be the same 
as if we simply measured the total oxygen abundance for 
the entire galaxy gas disk out to the radius R. For ref- 
erence, typical observed H I disks (e.g., Boomsma et al. 
2008 in NGC 6946) have radii of i? - 20 kpc. We consider 
a model spiral galaxy with a gradient of 12-|-log(0/II) of 
slope Fq/h units of dex per kiloparsec; Fq/h is and is 
typically of order — 0.05 dexkpc"^ (Zaritsky et al. 1994; 
van Zee et al. 1998; Kennicutt et al. 2003). We wilU also 
assume that our model galaxy has a gas surface den- 
sity S that obeys a power-law relation with the radius 
such that S(i?) = T,o{R/Ro)°'. While a constant surface 
density (a = 0) is a reasonable assumption at large radii 
(Wong & Blitz 2002; Leroy et al. 2008), we will show that 
the calculated oxygen abundance dilution is not strongly 
dependent on the choice of this power-law slope. For no- 
tational simplicity, in this section only we will use O/H 
as shorthand for the more cumbersome 12 -I- log(0/H). 

Let (0/II)sDSS be the mean O/H within an SDSS spec- 
troscopic fiber of a spiral galaxy (i.e., within the central 
3" diameter). Formally, 



(O/H) 



SDSS 



/ S(i?)(0/H)i?di?d6i 



/ J:{R)RdRde 





(5) 



where i?3 is the radius in kiloparsecs of the galaxy at a 
radius of 1.5" (i.e., spanned by the 3" SDSS fiber; this 
corresponds to « 4.9 kpc at z = 0.2) and we are implic- 
itly assuming the scale-height of the galaxy is roughly 
constant with radius. For radii R > R3, the oxygen 
abundance is O/H = (0/H)o + Tq/hR, where (0/H)o is 
the oxygen abundance extrapolated to i? = using the 
adopted the abundance gradient. We can now expand 
Equation (5) to find 



(O/H) 



SDSS — 



2 + a 
2nRi+°' 



(0/H)a 



27r J i?"[(0/H)o + ro/H^]i?di? 

(6) 





2 + a 

3 + a 



Solving for the R = abundance, we find (0/H)o = 
(0/H)sDss - [(2 + a)/(3 + a)]To/HR3- The mean metal- 
licity (0/H)</j out to a radius R is therefore 



fT.{R){O/ll)RdRd0 



(0/H)<^ = 



= (0/H)c 



/ i:{R)RdRde 





O/H 



R. 



(7) 

Hence the change in metallicity A(0 /H) we would expect 
within an SDSS fiber diameter is A(0/H) = [{2+a)/{3 + 
a)]ro/HAi? = [(2 + a)/(3 + a)]To/H{R - Rs)- For an 
outer radius of i? = 20 kpc, an inner radius of i?3 = 5 kpc, 
an abundance gradient of Fq/h — — 0.05dexkpc~^, and 
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Fig. 11. — Residual in 12 + log{0/H) from Tremonti et al. (2004) fit to the mass-metaUicity relation v. specific star formation rate 
(SSFR), both within the SDSS 3" fiber (left) and aperture-corrected (right). In each panel, the blue points are the high-mass low-metallicity 
outliers in our sample while the large grey points are the metal-rich dwarf galaxies from Paper I. (The "§" symbols denote the only two 
clearly spiral galaxies in the low-metallicity high- mass sample.) The grey cloud of points are the galaxies from the main Tremonti et al. 
(2004) SDSS sample described in §2, and the central thick black histogram is the median residual 12 + log(0/H) in bins of SSFR. The thin 
black histograms at the top and bottom of each panel are the distribution of SSFR of the parent sample, while the blue histograms at the 
bottom of each panel is the SSFR distribution of the low-metallicity outliers and the grey histogram at the top of each panel is the SSFR 
distribution of the high-metallicity outliers. 



a flat radial gas-surface density distribution of a = 0, 
this leads to an oxygen abundance dilution of A(0/H) = 
— 0.5dex, which is in the range we observe for these low- 
metallicity outliers (see, e.g., Fig. 11). 

We note that there are several caveats with this pic- 
ture. First, we have assumed that the observed oxygen 
abundance corresponds to a mass-weighted abundance, 
whereas because we measure the metallicity via emission 
line spectra, it is actually weighted by the luminosities of 
the H II regions. Secondly, it is unclear if the timescales 
involved pose any problem. Simulations show mergers 
do induce gas inflow to the centers of galaxies (Cox et al. 
2006), but most of this inflow occurs during the first close 
encounter between the two galaxies. While a few of the 
galaxies in our sample do have nearby neighbors, most do 
not: SDSS detects only a single, if disturbed, galaxy. It is 
possible that these are not yet dynamically relaxed; the 
two galaxies can pass close enough to one another that 
they are seen by SDSS as one object. If a merger-induced 
starburst lasts ~ 1 Gyr (Cox et al. 2006), then ^ 1% of 
all (properly defined) starburst mergers should display 



lower-than-expected metallicities. The relevant dynami- 
cal timescale for gas inflow depends on the typical radius 
from which the gas originates. Because the dynamical 
time idynam ~ -R/vcsc, gas at larger radii takes longer to 
reach and help dilute the central metallicity. 

Several studies have shown that galaxies in close pairs 
tend to have lower metallicies than expected. Kew- 
ley et al. (2006a) found that galaxies separated by < 
20h~^ kpc have systematically lower oxygen abundances 
by - 0.2 dex. On the other hand, Ellison et al. (2008b) 
found an offset of ^ 0.05 dex in 12 -I- log(0/H) at flxed 
mass galaxies in pairs with similar separations; the in- 
vestigations of Michel-Dansac et al. (2008) of galaxies in 
close pairs suggest that the relative log(0/H) offset is 
greater for the lower-mass galaxy in the pair than for 
the higher-mass galaxy. Kewley et al. also found rela- 
tively high speciflc star formation rates for the lower- 
metallicity member in their galaxy pairs; they attribute 
both the strong starburst activity and the lower oxygen 
abundance to fresh gas inflowing to the galaxies' central 
regions due to interaction with the closely neighboring 
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galaxy. 

In a similar vein, Lcc ct al. (2004) found that mor- 
phologically disturbed galaxies have systematically lower 
oxygen abundances (or higher luminosities) than galax- 
ies of similar brightness (or metallicity) . Likewise, lumi- 
nous infrared galaxies (LIRGs), whose high specific star 
formation rates and highly disturbed morphologies are 
believed to be due to a recent major merger, are ob- 
served to have oxygen abundances that are ~ 0.4 dex 
lower than other galaxies of comparable masses (Rupke 
et al. 2008). The pronounced blue colors of our outliers 
(see Fig. 9) implies that while the star formation rates 
for many of these galaxies are consistent with those of 
LIRGs, these outliers are clearly not as dusty as their 
infrared-luminous cousins; it is unclear the extent to 
which these classes of galaxies are related, and why, if 
LIRGs do typically have such relatively low metallicities, 
no galaxies in our sample are extremely red. 

4. CONCLUCIONS AND IMPLICATIONS 

We have identified a sample of 42 low-metallicity high- 
mass outliers from the mass-metallicity relation. As a 
population, these galaxies have disturbed morphologies 
and high specific star formation rates, implying that they 
are undergoing a merger-induced starburst. We propose 
that their observed low oxygen abundances are due the 
tidal interaction inducing large-scale gas inflow which 
subsequently dilutes the central interstellar medium in 
these galaxies. While there have been several observa- 
tional and theoretical suggestions that interacting galax- 
ies will result in a decreased observed gas-phase metal- 
licity, this is the first work that has shown that the vast 
majority of severe low-metallicity outliers from the mass 
metallicity relation are morphologically disturbed. Fi- 
nally, we note that the cuts outlined in §2.1 provide 
an effective means of using only colors and spectra to 
identify a rather pure (though not necessarily complete) 
sample of tidally disturbed galaxies. 

One striking implication of these results is that, while it 
is safe to assume that the metallicities for populations of 
galaxies will fall within a particular range of values given 
their luminosities at redshifts for which the luminosity- 
metallicity relation has been measured, one should not 
assign metallicities to individual galaxies based solely on 
their luminosities. For example, while studies of the 
host galaxies of long gamma-ray bursts (GRBs) suggest 
that GRBs are only found in low-metallicity environ- 
ments (Stanek et al. 2006; Kewley et al. 2007), several 
recent GRBs have been associated with very luminous 
hosts, such as GRB 070306 and its Mb ~ -22.3 host 
galaxy (Jaunsen et al. 2008). As we have shown in this 
paper and Paper I, assuming an oxygen abundance for 
a galaxy from its luminosity alone can result in mis- 
estimating 12H-log(0/H) in luminous galaxies by a whole 
dex, and thus one should not assume that, e.g., the host 
of GRB 070306 necessarily has a high metallicity. For 
cases in which log(0/II) has been measured via emission 
line spectra for GRB hosts, it is generally low, even if 
the galaxy is luminous. For example, the host galaxy of 
GRB 031203 has an absolute i?-band magnitude simi- 
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lar to that of the Milky Way, yet Margutti et al. (2007) 
find it to have a low metaUicity (12 -t- log[0/II] = 8.12) 
as well as a high star formation rate (~ 13Moyr~^). 
Prochaska et al. (2004) interpret this offset in metallic- 
ity as a sign that GRB 031203 went off in a "very young 
star-forming region," which is in line with our interpreta- 
tion in § 3.2. Futhermore, the brighter GRB hosts stud- 
ied by Fruchter et al. (2006) are morphologically very 
similar to our massive low-metallicity galaxies, and our 
results strongly suggest that the oxygen abundances in- 
ferred from luminosity alone are uncertain by as much 
as 1 dex. Specifically, if a luminous galaxy is morpho- 
logically disturbed, has a high specific star formation 
rate, and is extremely blue, then it should not be as- 
sumed to lie within the luminosity-metallicity locus of 
star-forming galaxies. 
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